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ABSTRACT 

In the last few years various methods of identifying struc- 
tural dynamics models from modal testing data have ap- 
peared. This paper presents a comparison of four of these al- 
gorithms: the Eigensystem Realization Algorithm (ERA), the 
modified version ERA/DC where DC indicates that it makes 
use of data correlations, the Q-Markov Cover algorithm, and 
an algorithm due to Moonen, DeMoor, Vandenberghe and 
Vandewalle. The comparison is made using a five mode 
computer model of the 20 meter Mini-Mast truss structure at 
NASA Langley Research Center, and various noise levels are 
superimposed to produce simulated data. The results show 
that for the example considered ERA/DC generally gives 
the best results; that ERA/DC is always at least as good as 
ERA which is shown to be a special case of ERA/DC; that 
Q-Markov requires the use of significantly more data than 
ERA/DC to produce comparable results; and that in some 
situations Q-Markov cannot produce comparable results. 

INTRODUCTION 

In recent years several time domain system identification 
techniques from the control theory community have become 
useful in structural modal testing, and still others have 
appeared which might offer advantages for modal testing. 
The purpose of this paper is to perform a comparison of 
four of these closely related identification methods on a 
typical structural dynamics identification problem. The aim 
is to make an initial determination of what the advantages 
and disadvantages of each method are, and under various 
different conditions, determine which method is best to use. 

In recent years, the Eigensystem Realization Algorithm 1 - 2 


1 Research Associate, Department of Mechanical 
Engineering, Old Dominion University, Member AIAA. 

2 Principal Scientist, Fellow, AIAA. 

3 National Research Council Senior Research Associate. 
Also, Professor of Mechanical Engineering, Columbia 
University, Fellow, AIAA. 


has seen considerable use on flexible spacecraft identifica- 
tion. This is one of the four closely related algorithms com- 
pared. An extension of this algorithm using data correla- 
tion (ERA/DC) is also studied. 3 The third approach is based 
on Q-Markov cover theory originally developed for model 
reduction. 4,6 It was generalized to identification using both 
white noise input tests and impulse input tests in Ref. 6, and 
then finally reformulated for tests using either type of input 
in Ref. 7. Although obtained from a different starting point, 
the resulting mathematical algorithm bears a family resem- 
blance to. ERA/DC. The fourth algorithm studied is that of 
Ref. 8 which allows direct use of general input tests. This 
algorithm will be refereed to as MDW, using the first letter 
of each author’s last name. A mathematical description of 
each algorithm is given in Appendix 1. 

THE TEST EXAMPLE 

One must use a mathematical or computer model rather 
than a real physical structural test article to compare the 
performance of the identification algorithms, since one must 
know what the correct identification result is. Here we pick a 
model obtained by finite element analysis of the Mini-Mast® 
structure illustrated in Fig. 1. This model is sufficiently 
complex to be a reasonable test of identification performance. 
The model for testing considers the first two bending modes 
as the lowest frequency modes with essentially identical 
frequencies, the first torsional mode, and the second two 
bending modes, again with essentially identical frequencies. 
Table 1 gives these frequencies and the associated damping 
expressed as the negative of the real part of the eigenvalues. 
The mathematical model is given in Appendix 2. 


Table 1 


Mode 

Damping 

Frequency 

No. 

Factor 

rad/sec 

HZ 

1 

0.09059 

5.032 

0.8008 

2 

0.09066 

5.036 

0.8014 

3 

0.32907 

27.420 

4.3640 

4 

0.38352 

38.351 

6.1038 

5 

0.38683 

38.682 

6.1565 
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The mathematical model considers two inputs and two 
outputs. The outputs are from nonorthogonal displacement 
sensors at the top of the structure, and the inputs are two 
torque wheels for the x and y axes, as shown in Fig. 1. 

THE APPROACH USED IN THE COMPARISON 

It is not a simple matter to determine a set of tests 
that constitute a reasonable comparison of the different 
identification approaches. Each approach has a somewhat 
different set of parameters that the user can adjust in order 
to get good performance. It would be of interest to optimize 
these parameters so that one compares in each case the 
best results each algorithm can produce. However, usually 
it is a matter of experience and intuition to decide how 
to make these parameter choices. In the case of ERA, 
methods based on perturbation theory have recently become 
available which help the ERA user know the accuracy of 
his identification results, and hence help him to adjust these 
parameters in order to improve the accuracy. 10-12 Similar 
results have been generated for ERA/DC and the Q-Markov 
cover approaches, and will be reported in the literature. No 
attempt is made here to optimize these parameter choices. 
The choices made are as follows. 

In Q-Markov, the measurement matrix is necessarily 
square, but in ERA and ERA/DC the matrix can be rectangu- 
lar, and in fact certain results from other sources 12 suggest 
that rectangular data matrices for ERA give bcuer results. 
Nevertheless, for ease of comparison, we restrict ERA and 
ERA/DC to square data matrices and recognize that their 
identification accuracy may be degraded. Also, in ERA/DC 
there is considerable freedom in the choice of the correlation 
matrix U(q), but here we restrict ourselves to the special case 
in which U(q) — H(q)H T (0) (see Appendix 1). During the 
comparison the size of the singular value decomposition ma- 
trix for each algorithm is kept the same (at 120x120), with 
the idea t hat the sin gular value decomposition in each al- 
gorithm is the most time consuming step. However, this 
makes the comparison of ERA with ERA/DC and with Q- 
Markov somewhat biased against ERA. Both ERA/DC and 
Q-Markov involve a parameter, s and d respectively, which 
determines the upper limit on a possibly large summation of 
products that does not appear in ERA. Hence, with the same 
computation time, ERA could have used a larger data matrix 
which could produce improved identification results. Some 
computer runs were made to study this. Appendix 1 gives 
a summary of each of the algorithms considered. Note that 
equal values of the parameters s and d in ERA/DC and Q- 
Markov imply essentially equivalent amounts of computation 
in setting up the respective data matrices. 

The lest example considered is a multiple-inpul and 
multiple-output one. As a result, it is of interest to have a 
scalar measure of the accuracy of the identification results 


to aid the comparisons of the various approaches. We 
accomplish this by considering the maximum singular value 
of the transfer function matrix of the identification error as a 
function of frequency. In other words, we take the maximum 
singular value at each frequency, of the true transfer function 
matrix of the computer test example minus the corresponding 
transfer function matrix of the identified system model. The 
symbol 9 is used to indicate maximum singular value. This 
scalar measure gives a simple way to express the combined 
error in natural frequencies, damping and mode shapes of 
the identified model. Note that this error criterion is used 
to design the model uncertainty weighting matrix for Hqo 
control design in Ref. 13. 

RESULTS 

Noise FYee Results 

Figure 2 gives identification results for, ERA, ERA/DC, 
and Q-Markov when the pulse response data is noise free. 
Figure 2a gives 9 of the true system transfer function matrix 
for comparison, and the 9 of the error transfer function 
matrix for ERA, and for ERA/DC with the parameters s 
chosen as 200, corresponding to 260 data points. The smaller 
the values of 9 of the error transfer function matrix curves, 
the better the identification results. Figure 2a shows that 
both ERA and ERA/DC give essentially perfect identification 
with the errors in the identification limited to the magnitude 
of numerical round-off. Figure 2b repeats 2a but for the Q- 
Markov algorithm. The similarity of the true and the d = 500 
curves in the second and third peaks means that these modes 
are not identified. 

We observe that the Q-Markov algorithm requires the use 
of very large data sets (e.g. d significantly greater than 2000) 
before it can approach the accuracy of ERA(with s equal 
59) and ERA/DC(with s equal 200) on noise free data. The 
same obs ervation can be seen in a different manner in Fig. 
2c. Only the first ten singular values should be nonzero, but 
the remaining singular values for the Q-Markov algorithm 
do not reduce to the level of numerical round-off until d is 
very large. 

Interpretation of Plots of Maximum Singular 
Value of the Error Transfer Function Matrix 

We are using the maximum singular value of the model 
error transfer function matrix, as a function of frequency, 
as our scalar measure of the model error. It is of interest 
to interpret these plots in terms of errors of frequency and 
damping of the modes in the model. In Figure 3a we show 
the maximum singular value of the error tran sfer f u n ction 
matrix when the model and the true system are identical 
except that the damping (i.e. the negative of the real part 
of the eigenvalues) of the first two modes is decreased by 
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0.1%, 1%, and 10%. The same computation was repeated 
using frequency error instead of damping error. When the 
frequency errors are 10 -s %, 10~ 2 %, and 10“ l % we obtain 
Fig. 3b which looks nearly identical to that of Fig 3a but with 
slightly lower peaks. Hence, the o(w) peak is approximately 
100 times more sensitive to frequency than it is to damping 
error, and the signatures of such errors are very similar. 

The Relationship Between ERA and ERA/DC 

In Appendix 3A we prove that ERA is a special case of 
ERA/DC in the following sense. Use the same Hankel matrix 
of pulse response data H( 0) in both algorithms. Pick the 
special case of U(0) = R{ 0) in equation (1.5) of Appendix 
IB for the singular value decomposition matrix in ERA/DC. 
Pick the same order k for the desired identified model from 
each algorithm. Then the pulse response history of the kth 
order ERA model is identical to the pulse response of the kth 
order ERA/DC model, for all choices of k, independent of 
the true order of the system. Hence the identified models are 
identical to within a linear transformation of state variables. 

Figure 2 a shows ERA/DC with s = 200. If s is set to 59 
in order to comply with the above assumptions in comparing 
ERA/DC to the ERA results, the a error curve in Fig. 2 a is 
similar to those shown, with differences again on the level of 
numerical round-off in the computations. This is consistent 
with the above theoretical result 

Note that unlike ERA, ERA/DC has the ability to in- 
crease the number of Markov parameters used in the iden- 
tification, without increasing the dimension of the singular 
value decomposition matrix, simply by increasing the num- 
ber of columns in H( 0). Since storage limitations constrain 
the size of the decomposition matrix, this extra freedom can 
be an important advantage of ERA/DC 

The Relationship Between ERA and Q-Markov 
Algorithms 

The relationship between ERA/DC and the Q-Markov 
cover algorithm is developed in Appendix 3B. ERA/DC uses 
a singular value decomposition of a matrix 11(0) in equation 
(1.6) to produce its realization, while Q-Markov uses singular 
value decomposition of the matrix D q in equation (1.11). 
Appendix 3B shows that D q can be produced from 11(0) 
by setting certain Markov parameters to zero in 11(0) (see 
equations 3.8 and 3.9 in Appendix 3B). 

As the parameter d in Q-Markov increases the number of 
Markov parameters set to zero becomes a smaller percentage 
of the total number of Markov parameters included in the 
identification. They also refer to later sample times which 
allows for more decay of the pulse response. 


The ERA/DC results for s = 59 and s — 200, which 
both have errors limited to the level of numerical round- 
off, suggest that ERA/DC gives accurate system realization 
from noise-free data for all r, s above the values needed 
to give row and column numbers of the Hankel matrix that 
are greater than the system order. The Q-Markov algorithm 
results suggest that it produces an accurate realization of 
the input-output relation only asymptotically as the amount 
of data used tends to infinity, i.e. d — ► oo. This is 
consistent with the analytical results in Appendix 3B; once 
d is sufficiently large that the Maikov parameters that are 
set to zero in the Q-Markov algorithm(equations 3.8 and 
3.9) become negligible, then Q-Markov produces results 
equivalent to ERA/DC with model errors reduced to the level 
of numerical round-off. This suggests that ERA/DC should 
give superior performance to Q-Markov Cover for all finite 
data sets. 

A Note on the Computational Effort of Each 
Algorithm 

The ERA, ERA/DC, and the Q-Markov Cover algorithms 
are compared by constraining the singular value decomposi- 
tion in each case to be of the same dimension, with the idea 
that the effort involved in this decomposition dominates the 
computational effort of the algorithm. Such a decomposi- 
tion is roughly a cubic function of the matrix dimension. 14 
ERA/DC and Q-Markov both require an additional multipli- 
cation of matrices to obtain the matrix which is decomposed. 
Such a multiplication is generally cubic in the matrix di- 
mension as well. However, the special structure of the data 
matrices allows one to generate a recursive procedure to sim- 
plify the computation. It produces a computation count that 
is quadratic in the matrix dimensions, so that the decompo- 
sition does in fact dominate. 

It is difficult to compare computational effort of the dif- 
ferent identification methods, because the details of the pro- 
gramming can be significant. Here we cite the CPU times 
used in running the different algorithms. All computations 
were carried out using MATLAB on a Micro VAX 3200 
workstation, running the high noise case discussed below. 
Table 2 gives the results for n = 120 (the number of rows 
of //,, or 11(d) in equation 3.9 of Appendix 3B, or 11(0)) 
in all cases, and for various values of s or d where s + 1 and 
cl -(- i arc the number of columns of Markov parameters in 
H(0) and 11(d) respectively. The actual number of columns 
is 2(s + 1) and 2(d+ 1), because the system has two inputs. 

Increasing s from 100 to 1500 in ERA/DC increased 
the CPU time by 24 seconds which represents a 26% 
increase in run time to include roughly 15 times as much 
data. The Q-Markov algorithm took slightly longer than 
ERA/DC for corresponding values of d and s. Also as d 
is increased from 100 to 1500 the CPU time increased by 
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40 seconds representing a 41% increase in run time. Note, 
however, that a much larger value of d is required in Q- 
Markov than of s in ERA/DC for comparable identification 
accuracy. Q-Markov with d - 1500 and ERA/DC with 
s = 500 give roughly comparable identification results, and 
Q-Markov takes 38% longer to run. We conclude that, at 
least for the programming implementations of the algorithm 
used, ERA/DC has a substantial advantage over Q-Markov 
in computer time needed. 


Table 2. Comparison of CPU Times for ERA, ERA/DC, 
and Q-Markov Cover ( n = 120 in all cases) 


s= 100 
92 


d- 100 
98 


ERA 


a= 59 

3 = 79 

3 = 99 


90 

120 

150 



ERADC 



3 = 300 

3 = 500 

s = 1000 

s = 1500 

96 

100 

110 

114 


Q-Markov 



a. 

II 

Co 

o 

o 

d= 500 

d = 1000 

d = 1500 

110 

117 

128 

138 


The ERA algorithm requires the singular value decompo- 
sition of a rectangular Hankel matrix 7/(0) rather than of a 
square matrix such as H( 0)^(0), one dimension of which 
is smaller. The increased dimension proportional to (s + 1) 
causes memory problems which prevented use of s = 300 
or above. In order to include additional data, one must con- 
vert to ERA/DC above this value. Since the CPU time in 
ERA is proportional to (s -f 1) we would expect a CPU 
time of 2252 seconds for s = 1500 if memory problems had 
not prevented such a computation. We conclude that ERA 
suffers from memory limitations, and that ERA/DC is to be 
preferred. Since ERA is a special case of ERA/DC in the 
prescribed sense given above, ERA/DC can always give as 
good results as ERA, and the extra parameters in ERA/DC 
can be optimized (e.g. by use of the methods of Refs. 10-12) 
to give better results. 


Low Noise Results 


is interesting to note in Figs. 4a and 4b that the drop in 
the singular values which is an indicator of the true system 
order, becomes worse (i.e. smaller) as either d or s is 
increased. ERA/DC is again better for corresponding values 
of d and $. This can be explained by examining the time 
constants for decay of each mode. For the highest frequency 
mode, four time constants is 10.34 sec which corresponds 
to a decay of the amplitude of the mode to 0.018 of its 
original value. When s = 800, the last data entry in the 
ERA/DC Hankel matrix corresponds to 25.8 sec, well after 
the highest frequency mode has disappeared into the noise 
level. Including more data points adds only noise and no 
signal as far as identification of the high frequency modes are 
considered. This added noise is reflected in higher singular 
values associated with noise, i.e. singular values 11 and 
above. On the other hand, four time constants of the lowest 
frequency mode is 44.16 sec, and adding data above s = 800, 
i.e. after 25.8 sec, still helps the identification of this mode, 
as is reflected in the increase of the first singular value with 
s. Obviously, there are different optimal values for s for 
different modes, generally with smaller values of 3 for higher 
frequency modes. 

Figure 4e shows that both ERA/DC and Q-Markov using 
large data histories associated with s = 800 and d = 2000, 
somewhat outperform ERA with its much shorter 119 point 
data history. 

High Noise Results 

Figure 5 repeats Fig. 4 but with a noise level increased by 
a factor of 10. Figure 5b shows that the Q-Markov algorithm 
fails to exhibit the drop in singular values at the true system 
order of 10, for all choices of & On the other hand, ERA/DC 
exhibits a drop at the true system order for lower values of s. 
One can interpret this in terms of the damping in the system. 
The Q-Markov algorithm appears to inherently need more 
data points than ERA/DC, but as more data points are taken 
the decay of the true system modes causes the extra data for 
times after significant decay has occurred to have poor signal 
to noise ratio. This increases the level of the singular values 
associated with noise, and decreases the step change in the 
singular values at the true system order. 


Figure 4 introduces a low noise level in the data. The 
noise is modelled as additive, white, zero-mean and Gaussian 
with standard deviation equal to 8 x 10~ 9 . Note that the 
signal initial root mean square value is 1 x 10~ 6 . 

Figures 4b and 4d again show that the Q-Markov method 
needs a long data set, e.g. d = 1500, before all modes 
are clearly identified. Figures 4a and 4c show that the 
same level of identification error is achieved by ERA/DC 
with many fewer data points, e.g. for s = 400, a number 
too small for successful identification with Q-Markov. It 


The comparison of the methods in Fig. 5e shows that 
ERA/DC using s = 400, and Q-Markov using roughly four 
times as much data with d = 1500, produce comparable 
results. 

High Damping Results 

We consider the same high noise situation as in the 
previous section, but change the system by increasing the 
damping in each mode of the system. The new real parts of 
the eigenvalues are -0.6, -0.6, -1.5, -1.8 and -1.8 respectively. 
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This corresponds to decreasing the time constants for the 
envelope of decay of each mode from 11.04, 11.03, 3.04, 
2.61, 2.59 to 1.67, 1.67, 0.67, 0.56, 0.56. 

Because of the fast decay, lower values of s are appro- 
priate in ERA/DC so as not to include data after complete 
decay of the signal. ERA/DC succeeds in identifying all 
modes. Figure 6a is somewhat congested in the neighbor- 
hood of the second and third peaks, but one can observe that 
when s increases from 60 to 100 the model error for these 
peaks gets worse while the model error for the first peak, cor- 
responding to the first two modes that decay more slowly, 
remains essentially unchanged. Thus, s -- 100 is beyond the 
optimal value for the high frequency modes. 

Figure 6b shows that the Q-Markov cover algorithm is 
essentially unable to identify or gives very poor identification 
of the three highest frequency modes. Unlike ERA/DC, the 
algorithm needs a d large enough that the Markov parameters 
that are set to zero either have sufficiendy small values due 
to the decay of the signal or are insignificant in number 
relative to the total number of other Markov parameters. 
Neither situation occurs before d becomes so large that 
increasing d simply adds noise without adding signal. The 
best compromise appears to be d = 200, but this is still 
insufficient for good identification of these modes. In this 
high noise and fast decay situation, ERA/DC can successfully 
identify all modes and Q-Markov Cover cannot 

Q-Markov Cover Results Using Random Input 
Response 

The Q-Markov Cover algorithm has the potential ad- 
vantage that it can use random input response histories di- 
rectly, rather than using data from such an input to obtain 
the Maikov parameters and then using these to perform the 
identification as described above. Figure 7 shows the results 
for numerical tests with noise-free data. The method iden- 
tifies the first two modes with comparatively poor accuracy, 
and totally fails to identify the remaining three modes. Go- 
ing above d = 1500 produces very modest improvement. 
Averaging a set of runs to obtain Maikov parameters and 
covariance parameters can produce additional modest im- 
provement. At least for the example considered, this poten- 
tial advantage of Q-Markov does not appear to be realizable 
in practice. 

MDW Results 

In Fig. 8 we tum to the MDW algorithm which is an- 
other algorithm that can use random input responses without 
the need to obtain Markov parameters first. In fact, it can 
use any sufficiently rich inputs. Like ERA, this algorithm 
uses a data matrix that need not be square. The originator of 
the algorithm suggests use of a rectangular matrix. Also, the 


matrix involved, which undergoes a singular value decom- 
position, contains not only the measurement data but also 
the input function. This means that one has the freedom to 
use general inputs as long as they are sufficiently rich, rather 
than obtaining impulse responses. This benefit is obtained at 
the expense of, in our case, doubling the size of the singular 
value decomposition matrix for the same amount of data. 

In place of the 120x120 matrices used in all previous 
computations the matrix size for H in these runs is chosen 
as 80x160, i.e. a rectangular matrix with a roughly compa- 
rable computation time. The inputs were taken as random, 
producing displacement sensor outputs shown in Fig. 8a. 
Figure 8b gives the singular values of H. The first 40 are 
associated with the random input functions, and the next 10 
are associated with the system dynamics. Experience indi- 
cates the input signals need to be very rich to obtain good 
results. Noise free results are given as well as results with 
three noise levels. The two higher noise levels of 5 x 10 - ^ 
and 5 x 10 -8 correspond to the high and low noise level of 
the previous computations, and the third level is associated 
with a decrease by another factor of 10. 

Like ERA and ERA/DC the identification error is on the 
order of numerical round-off error in the no-noise case. Fig. 
8c. Figure 8d indicates somewhat worse error behavior than 
the other algorithms when noise is present. 

CONCLUSIONS 

Some of the conclusions suggested by the numerical test 
results are as follows: 

1. Of the four algorithms tested ERA/DC gives the best 
results. 

2. ERA is shown to be a special case of ERA/DC in the 
sense that in the no noise case and using the same 
Hankel matrix for each method, the models produced by 
the two algorithms give identical input/output relations 
for all model orders chosen (i.e. for all choices of the 
number of nonzero singular values retained). The state 
space realizations for the models can be different, and 
there can be different numerical round-off properties, 
but otherwise, they both give the same input-output 
relationship. 

3. They both produce models with identical input-output 
relations when each uses the same Hankel matrix, but 
ERA/DC has a computational advantage since it produces 
a smaller singular value decomposition for rectangular 
Hankel matrices. In practice, memory limitations asso- 
ciated with the singular value decomposition dictate the 
maximum size of the rectangular Hankel matrix in ERA, 
based on both the shorter and longer dimensions of the 
matrix. In ERA/DC the singular value decomposition is 
of a square matrix whose dimension is the shorter dimen- 
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sion of the Hankcl matrix. This allows ERA/DC to give 
significantly better results in practice. 

4. The Q-Markov Cover algorithm can be obtained from 
ERA/DC by setting a triangular set of Markov parameters 
to zero in the data matrix, and then using a different 
realization. As a result, Q-Markov is the only algorithm 
tested here that needs large sets of data to obtain accurate 
identification from noise-free data. The setting of these 
Maikov parameters to zero can be insignificant if they 
apply to times after all transients have decayed, or if 
enough Maikov parameters are included that the set 
becomes insignificant. Hence, Q-Markov required much 
longer sets of data to be used (compare s and d values 
in the figures) in order to obtain results comparable 
to those of ERA/DC. Also, situations can arise where 
the speed of decay of the transients prevents Q-Markov 
from identifying a mode that ERA/DC can successfully 
identify. 

5. The Q-Markov Cover algorithm has the potential advan- 
tage that in one of its forms it can accept random input 
response data directly. ERA and ERA/DC, as well as 
another form of Q-Markov, require computing pulse re- 
sponse histories from the data associated with the actual 
inputs, before the algorithms can be applied. Our exam- 
ple run with Q-Maikov using random input response data 
directly, gave significantly worse results than the other 
algorithms. This suggests that the potential advantage is 
not actually an advantage in practice. 

6. The MDVV algorithm can directly use data from any 
sufficiently rich input signal, but care must be exercised 
to insure a well conditioned full rank input signal. This 
advantage is counterbalanced by the disadvantage of 
requiring a singular value decomposition of a larger 
matrix than all of the above algorithms when the same 
amount of data is used. Results appear to be inferior to 
the other methods. 

We reiterate that the conclusions presented here are based 
on testing only one structural dynamics example, and that 
the comparisons are made without fully optimizing all the 
parameters of any of the methods. In practice, the variance 
and bias computation methods developed using perturbation 
methods (Refs. 10-12) can be used with experimental data 
to help the user of ERA, ERA/DC, or Q-Markov Cover to 
optimize the parameters of these methods. 
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APPENDIX 1 

MATHEMATICAL DESCRIPTION OF THE 
IDENTIFICATION ALGORITHMS 


A. ERA Algorithm 


1. Form the Hankel matrix 


H(j~ 1 ) = 



J+2 


a 


Y i +*+ 1 


( 1 . 1 ) 


Y j+r + 1 Y j+r+t / 


where is the ith Markov parameter of the pulse 
response history. 

2. Form the singular value decomposition of 77( 0) 


77(0) = U'EV 7 ' ( 1 . 2 ) 


3. The system matrices of the ERA identified kth order 
discrete time model are given as 

At = E 77(1)HE " 1/2 


5 fc = Ej / 2 Vf£^ (1.3) 

= 7f U^J 2 


B. ERA/DC Algorithm 

1 . Form the Hankel matrix as in equation (1.1) 

2. Compute Rfq) 


tyq) = Hiq^i 0) (1.4) 


3. Form block correlation matrix 


U(q)= 

Biq) R(q + 7 ) 
77(9 + 7 ) RlQ + 27 ) 

\R(g + Q7) 


R(q + P 7 ) \ 

77(9+(o + /3)7)/ 
(1.5) 


4. Form the singular value decomposition of Z/(0) 


7/(0) = UZV T ( 1 . 6 ) 


5. The ERA/DC identified kth order discrete time model is 
A* = E- 1/2 t^7/(l)F*E; l/2 

/ ^( 0 ) \ 

B k = e; 1/ V (1.7) 

\ 77 ( 07 )/ 

c k = £?u k e;/ j 


C. Q-Markov Cover Algorithm 
1. Calculate X, 


d 


X,- = Y;(c7,i) = 

k = 0 


(1.8) 


The matrix E* is the upper left hand kxk partition of 
E containing the k largest singular values along the 
diagonal. Matrices U k and V k are obtained from U and 
V by retaining only the k columns, i.e. singular vectors, 
associated with the singular values retained in E*. Matrix 
Ei is a matrix of appropriate dimension containing l 
columns, that is all zero except that the top 7x7 partition 
is an identity matrix. Em is defined analogously. 


here Yi is the pulse response at time i. 
2 . Form Rq, H q and Mq 

/ Xo XJ ... 



U’-i *0 / 
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Y 0 0 ... 0 \ / 

H q = Yl y ° •” ; , M q = 

W-x YoJ V 

3. Compute D q 

D q = 1^ - H q Hj 

4. Form the singular value decomposition of D q 


D q = UEV T 


( 1 . 9 ) 

( 1 . 10 ) 


( 1 . 11 ) 


5. Form P k and partition it 


Pk = u k z 1 k / \ P k = } , Oi e ir* 


( 1 . 12 ) 


6. Form P and P 


( Oq \ / Oi 

°'lp=h 


( 1 . 13 ) 


7. The Q-Markov identified kth order discrete model sys- 
tem matrices are 

A k = P + P , B k = Oq, C k = P + M q ( 1 . 14 ) 

where + indicates the Moore-Penrose pseudoinverse 

D. MDW 

l. Form Hy, and II 


u k+i Uk+i+l 
Vk+i Vk+i + 1 

u k+i + 1 u t+i+2 
2/*+«'+l Vk+i+2 

u k+2i-l «k+2i 
,Vk+2i-l Vk+2i 


u k+i+j - 1 
Vk+i+j-l 
u k+i+j 
Vk+i+j 

u k+j+2i-2 

yk+j+ 2%-2 , 


(I) 


( 1 . 15 ) 


where u k and y k denote l dimensional input vector and 
m dimensional output vector at time k. 

2. Form the singular value decomposition of H 

H=lKVr = (Vn (U«) 

here the sizes of the matrices are 

size(Un) = (mi + li) X (2 mi + n) 
size(Ui 2 ) — (mi + li) x (2li — n) 
size(U 2 \) = (mi -f li) x (2 mi 4- n) 
size(U 22 ) = (mi -b li) X (2/t — n) 
size(Sn) = (2 mi -f n) X (2mi + n) 

3. Form the singular value decomposition of UtfUnSn 

lf[,U u S u = [U, o)(vS-) (1-17) 

where U q is a (2/i - n) x n matrix. 

4. The identified nth order model 

(U?U&U(m + l+l : (i+ l)(m+/),:)S\ 
\U(mi + li+m+l:(m+l)(i + l),:)S J 


A B 
C D 


UfUTU{l:{m+l)i,:)S \ 
U(ttu + li -I- 1 : mi + li + m, :)S J 

( 1 . 18 ) 


APPENDIX 2 

MINI-MAST MODEL 


«k 

«k+l • 


Vk 

Vk+l ■ 

Vk+j-l 

Ufc+l 

u k+ 2 ■ 

Uk+j 

Vk+i 

Vk+2 ■ 

Vk+i 

u k+i - 1 

Uk+i ■ 

■ ■ u k+j+i- 

Vk+i-l 

Vk+i • 

■ ■ Vk+j+i- 


x = Ax + Bu 
y-Cx 
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A - 


Columns 1 through 5 


- 0.0892 

- 0.0001 

- 5.0354 

0.0000 

0.0009 

0.0001 

- 0.0891 

- 0.0001 

5.0316 

0.0131 

5.0354 

0.0002 

- 0.0921 

- 0.0001 

- 0.0014 

0.0000 

- 5.0316 

0.0001 

- 0.0920 

- 0.0154 

- 0.0032 

- 0.0021 

0.0035 

0.0074 

- 0.3251 

0.0040 

- 0.0076 

- 0.0040 

0.0027 

27.4201 

0.0247 

- 0.0954 

- 0.0269 

- 0.1040 

0.0015 

- 0.0956 

- 0.0251 

0.1042 

- 0.0275 

0.0038 

0.0266 

- 0.1015 

- 0.0262 

- 0.0997 

- 0.0033 

- 0.1020 

- 0.0263 

0.1005 

- 0.0257 

0.0149 

Columns 

6 through 

10 



- 0.0015 

0.0055 

- 0.0200 

- 0.0099 

0.0374 

0.0153 

- 0.0211 

- 0.0061 

0.0386 

0.0108 

0.0012 

0.0101 

- 0.0377 

- 0.0055 

0.0205 

- 0.0141 

0.0388 

0.0108 

- 0.0218 

- 0.0065 

- 27.4201 

- 0.0094 

0.0177 

0.0111 

- 0.0252 

- 0.3330 

- 0.0099 

0.0224 

0.0112 

- 0.0212 

0.0028 

- 0.3762 

0.5972 

38.3637 

- 0.1010 

- 0.0132 

- 0.5956 

- 0.3790 

- 0.0466 

38.6603 

- 0.0028 

- 38.3637 

0.0464 

- 0.3912 

- 0.5969 

- 0.0027 

0.1011 

- 38.6603 

0.5986 

- 0.3943 


B - C - 


0.00235 

- 0.00200 

0.01119 

- 0.00911 

- 0.00201 

- 0.00236 

0.00402 

0.00762 

- 0.00235 

0.00200 

0.01122 

- 0.00914 

- 0.00201 

- 0.00236 

- 0.00402 

- 0.00764 

- 0.00011 

- 0.00025 

- 0.00917 

- 0.00931 

0.00011 

0.00025 

- 0.00918 

- 0.00933 

0.00167 

0.00095 

- 0.00043 

- 0.00243 

- 0.00091 

0.00155 

- 0.00245 

0.00196 

0.00163 

0.00092 

0.00047 

0.00242 

- 0.00089 

0.00151 

0.00239 

- 0.00199 


APPENDIX 3 

COMPARISON OF IDENTIFICATION 
ALGORITHM 

A. Comparison of ERA and ERA/DC 

Wc can show that if one uses the same Hankel matrix 
//( 0) in ERA and in the special case of ERA/DC with 
M(0) = tf(O), then the resulting realizations, for any chosen 
system order, have identical pulse response histories. In 
other words ERA is a special case of ERA/DC in the 
sense that the realizations produced have the same Markov 
parameters, but they can have different choices for state 
variables. 

The (j + 1) Markov parameter of the kth order model 
obtained by ERA in equation (1.3) is 


= £?U k X k V?UlH(l)V k yV?Em (3.1) 

Construct the (j + 1) Markov parameter of the kth order 
model from ERA/DC using the W(0) = R{ 0) = H{ 0)^(0) 
special case. Use V , E, and V for the singular value 
decomposition of H( 0) as in (3.1), and use V, E, V for the 
singular value decomposition of equation (1.6) for ERA/DC. 
Then these decompositions are related by 

0=U, E = EE t , V=U (3.2) 

Substituting these decompositions into (1.7), produces the 
realization 

A k = x; 1 ulW)U k x; 1 

B k = V- k 1 UlimEm (3.3) 

C k = E?U k Efc 

whose ( j + 1) Maikov parameter is 

Yj+i =£f u k v k {Y.- k l if£ H(i)VE T u T u k i:; 1 y 

X^UlUTA^Em (3.4) 

and after noting that 

VTflFU&i 1 = V k (3.5) 

One obtains equation (3.1), which completes the proof. 

B. Comparison of ZY(0) in ERA/DC and Dq in 
Q- Markov 

ERA/DC creates its realization from a singular value 
decomposition of W( 0). We pick the special case ZY( 0) = 
0) = H(0)H T (0) in this paper. The i,jth partition of 
M(0) can be written as 

im)ij = £ *'.j = l»2,...,r+l (3.6) 

Jb =0 

The Q-Markov algorithm obtains its realization from a sin- 
gular value decomposition of the symmetric matrix Dq. The 
ijth partition for i > j can be written as 

< i+i i~ l __ 

mu = YYi-j+kYl ~ Y/i-i-tf-i -h 
k = o t=o 

d + l 

= y3^*-j+k^T - y 'Yj- j+ k Y k 

k = 0 k=0 


Y j+i = C k A> k B k 
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d+1 

= E 1 

d+i-y 

= ^ Yi + ifYj+ k i y j = 1,2,. . . ,q (3.7) 

k=0 


Comparison of these two formulas shows that if we pick 
s — d, and q = r + 1, then the first row partition of D q , 
i.e. for j = 1, is identical to the first row partition of 14(0). 
However, for larger j there are j — 1 terms missing in the 
summation for [D 9 ],y that are present in [£/(0)],y. 

An alternative way to view the difference is to write the 


following matrix expressions 

1/(0) = H(0)H t (0 ) 

[Dq)ij = (H (0)HJ )ij for i>j 


where 

( Yl 

Y, . 

^+1- 


Y 2 

Yi . 

.. Y d+ 2 _ 

: ; 

; 


U 1 - 

•• Y d+l 


0 


(3.8) 


Y d 

Y d+ i 


Y d+ i\ 
0 


.. 0 / 

(3.9) 

and 7T(0) is given in equation (1.1). Thus, the 1/(0) differs 
from D q by setting the lower right triangular set of Markov 
parameters in the H( 0) to zero in order to produce H d . 



Tip Plate 


y Torque Wheel 


-X 


x Torque Wheel 


y 



Figure 1 Mini-Mast test article showing the x and y torque wheel locations and 
the A, B displacement sensor locations and their axes on the bay 18 tip plate. 
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Frequency (Rad/Sec) 


Frequency (Ratl/Sec) 


Figure 2 Noise free identification results: 

(a) Maximum singular value of true transfer function matrix vs. frequency, 
and of the model error transfer function matrices for ERA and ERA/DC 
(s = 200) vs. frequency. 

(b) Maximum singular value of true transfer function matrix vs. frequency, 
and of the Q-Markov model error transfer function matrices vs. frequency for 
various d . 

(c) Singular values of measurement matrix for Q-Markov identification for var- 
ious d. 



Frequency (Rad/Sec) Frequency (Rad/Sec) 

Figure 3 Maximum singular value of the model error transfer funetion matrix as a 
g function of frequency: (a)for various damping errors in the first two modes, 
(b)for various frequency errors in the first two modes. 


11 






3 = 


3 = G00 


8 


d = 1000 
d = 1500 
d = 2000 


3 = 

3 = 600 

5 = 800 


10 10 



Frequency (Rad/Sec) 


— </ = 500 

- - </= 1000 
... d = 1500 

</ = 2000 


— TVue : 

- - ERA I 
... Q- Markov . 

(d = 2000) | 
... _ ERA/DC : 
\ (5 = 800) ‘ 


Frequency (Rad/Sec) 


Frequency (Rad/Sec) 


Figure 4 Low noise identification results: 

(a) Singular values of measurement matrix for ERA/DC identification for var- 
ious 5, 

(b) Singular values of measurement matrix for Q-Markov identification for 
various d , 

(c) Maximum singular value of the ERA/DC model error vs. frequency for 
various s . 

(d) Maximum singular value of the Q-Markov model error vs. frequency for 

various d. ------ — — - 

(e) Comparison of maximum singular value of true transfer function matrix 
and of the model error transfer function matrices for ERA, ERA/DC (s = 
800), and Q-Markov (d = 2000). 











Frequency (Rad/Sec) Frequency (Rad/Sec) 

Figure 5 High noise identification results: 

(a) Singular values of measurement matrix for ERA/DC identification for var- 
ious 3. 

(b) Singular values of measurement matrix for Q-Markov identification for 
various d . 

(c) Maximum singular value of the ERA/DC model error vs. frequency for 
various s. 

(d) Maximum singular value of the Q-Markov model error vs. frequency for 
various d. 

(e) Comparison of maximum singular value of true transfer function matrix 
and of the model error transfer function matrices for ERA, ERA/ DC (s = 
400), and Q-Markov (d = 1500). 
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Frequency (Rad/Sec) Frequency (Rad/Sec) 


Figure 6 Maximum singular value of the model error transfer function matrix for fast 
decay system: (a)ERA/DC, (b)Q-Markov Cover. 



Frequency (Rad/Sec) Frequency (Rad/Sec) 


Figure 7 Maximum singular value of the model error transfer function matrix for Q- 
Markov Cover using random input response data directly. 
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xlO-s (a) (b) 





Frequency (Rad/Sec) Frequency (Rad/Sec) 

Figure 8 Identification results for the MDVV algorithm: 

(a) Time histories of the A and D simulated bay 18 displacement sensor out- 
puts for the chosen random iuput signals. 

(b) Singular values of the MDW input-output matrix with noise free data, 
and with three noise levels. 

(c) Maximum singular value of the true transfer function matrix and of the 
model error transfer function matrix (noise free) vs. frequency. 

(d) Maximum singular value of the true transfer function matrix and of the 
model error transfer function matrices with various noise levels vs. frequency. 
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